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Abstract 

Many plant populations have adapted to local soil conditions. However, the role of arbuscular mycorrhizal fungi is often 
overlooked in this context. Only a few studies have used reciprocal transplant experiments to study the relationships 
between soil conditions, mycorrhizal colonisation and plant growth. Furthermore, most of the studies were conducted 
under controlled greenhouse conditions. However, long-term field experiments can provide more realistic insights into this 
issue. We conducted a five-year field reciprocal transplant experiment to study the relationships between soil conditions, 
arbuscular mycorrhizal fungi and plant growth in the obligate mycotrophic herb Aster amellus. We conducted this study in 
two regions in the Czech Republic that differ significantly in their soil nutrient content, namely Czech Karst (region K) and 
Ceske Stredohori (region S). Plants that originated from region S had significantly higher mycorrhizal colonisation than 
plants from region K, indicating that the percentage of mycorrhizal colonisation has a genetic basis. We found no evidence 
of local adaptation in Aster amellus. Instead, plants from region S outperformed the plants from region K in both target 
regions. Similarly, plants from region S showed more mycorrhizal colonisation in all cases, which was likely driven by the 
lower nutrient content in the soil from that region. Thus, plant aboveground biomass and mycorrhizal colonisation 
exhibited corresponding differences between the two target regions and regions of origin. Higher mycorrhizal colonisation 
in the plants from region with lower soil nutrient content (region S) in both target regions indicates that mycorrhizal 
colonisation is an adaptive trait. However, lower aboveground biomass in the plants with lower mycorrhizal colonisation 
suggests that the plants from region K are in fact maladapted by their low inherent mycorrhizal colonization. We conclude 
that including mycorrhizal symbiosis in local adaptation studies may increase our understanding of the mechanisms by 
which plants adapt to their environment. 
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Introduction 

Plant populations often adapt to their local environments [1,2]. 
Local adaptations can be demonstrated by reciprocal transplant 
experiments showing that local genotypes perform better in their 
native environments than genotypes from other populations [3]. 
Theoretically, the magnitude of local adaptation should increase 
with increasing genetic variation within populations and with 
increasing environmental and phenotypic divergence between 
populations [2] . As a result, local adaptation is primarily expected 
at large spatial scales; however, it has also been found on small 
spatial scales (e.g., between microhabitats a few metres apart) [4] . 

Local adaptations in plants can often be explained by 
differences in the soU conditions of the studied localities [5]. 
While the effect of soil is most commonly explained as the effect of 
chemical or physical soil properties [6,7], it may also be related to 
the differences in soil biota [5]. Recentiy, several studies have 
suggested that plants may be locally adapted to arbuscular 
mycorrhizal fungi [8-11]. It has been shown that plants perform 
better when they grow in their local soil and when they are 
inoculated with their native arbuscular mycorrhizal fungi [8-10]. 



In contrast, another study of 64 plant species showed that the 
responses of plants to both native arbuscular mycorrhizal fungi 
and local soil can vary largely, from negative to positive effects 
[1 1]. One study indicated that plants can change their dependence 
on arbuscular mycorrhizal fungi depending on the actual nutrient 
content of the soil [8]. Other studies suggest that the ability of 
plants to change the intensity of mycorrhizal colonisation is limited 
[9,12]. 

The association of plants with arbuscular mycorrhizal fungi 
cannot be generally regarded as mutualism, but it can be found 
somewhere along the mutualism-parasitism continuum [13]. 
Mutualistic associations are most likely to arise in nutrient-limited 
environments, and parasitic associations are most likely in high- 
fertility environments [14]. Because the acquisition of soil 
resources can act as a strong selection pressure, it is likely that 
plants and their associated arbuscular mycorrhizal fiangi exert 
reciprocal selective forces on each other [12]. As a result, plants, 
arbuscular mycorrhizal fungi, and other soil organisms are selected 
to their ability to coexist within local communities under ambient 
soil conditions (called the co-adaptation model) [12,15]. In 
agreement with this theory, the effects of arbuscular mycorrhizal 
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fungi on the number of shoots and aboveground biomass of several 
vascular plants have been found to vary among sites with different 
soil fertilities [16]. 

Only a few studies have used reciprocal transplant experiments 
to study the local adaptations of plants to soil conditions mediated 
by mycorrhizal fungi. Furthermore, most of the studies were 
performed under controlled greenhouse conditions [8-10,17,18]. 
Nevertheless, controlled conditions cannot include all naturally co- 
occurring organisms such as pathogens and herbivores, and they 
also cannot include all of the abiotic factors of the sites, such as 
local microclimatic conditions. As a result, field reciprocal 
transplant experiments regarding local adaptation often produce 
different results from those in common gardens [7,19]. Field 
experiments are expected to provide more realistic insight into the 
issue of local adaptation than common garden experiments [19]. 

We studied the role of arbuscular mycorrhizal fungi in 
adaptations oi Aster amellus to local soil conditions. We conducted 
this study in two regions that differ significantly in their soil 
nutrient content [20,21]. A previous study indicated that this 
species is an obligate mycotrophic herb with high dependency 
upon mycorrhiza because all plants in the non-inoculated 
treatment died within a few months [20]. We expected to find 
local adaptation in A. amellus because its populations are 
genetically differentiated [21]. Furthermore, we found some 
evidence of local adaptation after two years in previous 
experiments with this species [21,22]. However, local adaptation 
was found only in the seedling stage and differed among 
populations [2 1] . We were interested in whether local adaptation 
could be detected after five years in the field and if this adaptation 
was influenced by the association of the plant with mycorrhizal 
fungi. We selected three sites in each of two regions with diploid 
populations of A. amellus. We collected seeds in the field and 
planted juvenile plants pre-grown from the seeds from each 
population into each population. We monitored plant growth in 
th(^ fiiJd for five years. After that, we determined the aboveground 
dry biomass and the percentage of mycorrhizal colonisation in the 
roots. We collected data in only four out of the six populations (two 
in each region). 

We asked the following questions: 1. Does the mycorrhizal 
colonisation of A. amellus plants differ among populations in the 
field? 2. Is there any evidence of local adaptation in A. amellus after 
five years in the field? 3. Do the differences in root colonisation 
between plants from different regions and populations match the 
patterns of plant growth in the field? 

Materials and Methods 

Ethics statements 

The administration of Czech Karst, a protected landscape area, 
permitted us to work with the endangered herb Aster amellus in the 
protected area and in the other locations. No specific permissions 

were required for field work in the Ceske Stredohori region 
because the sites were not privately owned or protected in any 
way. 

Study species and sites 

Aster amellus L. (Asteraceae) is an endangered, perennial plant 
growing in dry grasslands in Europe and Asia [23]. A. amellus 
occurs as a diploid and hexaploid cytotype in the Czech Republic 
[23]. The two cytotypes occur in close proximity, but they never 
form mixed populations in the Czech Republic [24] . We used only 
diploid populations for this study. 

We conducted this study in two regions in the Czech Republic 
that differ in their geological bedrock, vegetation composition and 



soil conditions. The first region, Ceske Stredohori (region S), is 
characterised by marl bedrock and vegetation belonging to the 
Bromion community [25]. The second region, Czech Karst (region 
K), is characterised by limestone bedrock and oak-hornbeam 
forests of the association Querco pubescenti-petraeae (thermophile oak 
wood) [26]. 

We selected three sites with populations of Aster amellus in each 
of the two regions (Table 1). The distance between the two regions 
was approximately 70 km and the distances between sites within 
regions ranged from 2.5 to 17 km. We measured the soil 
conditions in all sites. We collected five samples of the upper soil 
layer at each site. We transferred them to the laboratory, air-dried 
them, sieved them through 2 mm mesh and homogenised them. 
In the laboratory, we determined the actual pH and analysed the 
concentrations of Ca^^, Mg^"'", following Moore and Chapman 
[27] and evaluated the available P using the photometric method 
[28], the total N and total C content following Ehrenberger and 
Gorbach [29] and the organic C and carbonate content using the 
volumetric method according to ISO-Standard (10693). We tested 
for differences in soil properties among sites and regions with a 
hierarchical analysis of variance (site nested within region) 
(Table 1). The nutrient content of the soil was generally low in 
all sites, but the sites in region K had significantiy more potassium, 
nitrogen and organic carbon contents but significantiy lower 
carbonate content than the sites in region S (Table 1). 

Five arbuscular mycorrhizal fungi species were found in the soil 
from the studied sites (Glomus mosseae, G. intraradices, G. microag- 
gregatum, G. etunicatum, and G. constrictum; H. Pankova and K. Krak, 
unpublished). Previous garden experiments and field observations 
indicated that the percentage of mycorrhizal colonisation differs 
between populations of ^. amellus, has a genetic basis and is likely 
driven by the nutrient content in the soil of the natural populations 
[10,20]. Mycorrhizal colonisation was found to be significantiy 
lower in region K than in region S [20]. 

Reciprocal transplant experiment 

We collected seeds of A. amellus from each of the six sites in 
September 2003. In February' 2004, we sowed c. 300 seeds from 
each site in plastic trays with garden substrate and placed them in 
the greenhouse at 10 °C with natural light conditions. In April 
2004, we transplanted the seedlings into multipot trays with pots of 
3x3 cm and kept them in a common garden. In October 2004, we 
transplanted small juvenile plants into the field. We conducted a 
full reciprocal transplant experiment, so that plants from each site 
were planted at each site. Therefore, we distinguish the target site 
(where the plants were planted experimentally) and the site of 
origin (where the plants had come from). At each target site, we 
used a randomised block design and planted 150 juvenile plants (6 
sites of origin x25 replicates) into five randomly placed rows. The 
distance between the two nearest plants in a row was 10 cm. We 
marked each plant with a nail and a metal label with a number for 
later recovery with a metal detector. We watered the plants 
immediately after planting. 

Because this perennial plant grows very slowly in the field, we 
harvested the experiment after five years, in October 2009. We 
recorded the survival percentage, the number of leaves and the 
number of flowers in the field. Because few were flowering (< 1 %), 
we could not analyse plant reproduction data in this study. We 
extracted each plant with large soil monoliths (approximately 
15 cm in diameter and 10 cm deep) and harvested both the 
aboveground biomass and the roots for laboratory analyses. In the 
laboratory, we weighed the aboveground biomass after drying it to 
a constant weight at 70 °C for 72 hours. We carefully washed the 
roots from the soU, cut them into segments and stained them with 
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0.05% Trypan blue in lactoglycerol [30]. We quantified the 
mycorrhizal colonisation in the roots using the modified segment 
method [3 1] under a compound microscope at 200 x magnifica- 
tion. 

Because of soil erosion at two sites and the subsequent damage 
of large parts of the experiment, we collected data from all six sites 
of origin only in four target sites (two in each region; Kl, K3, SI, 
S2). In those four sites, we planted 600 plants in total. However, 53 
'm plants were destroyed by animals at the beginning of the 

0 experiment. Therefore, we calculated the percent survival out of 
^ 547 plants. The mean survival percentage was 5 1 % (Table 2), 

which is quite high after five years in the field. Thus, the data on 
" plant size (number of leaves and aboveground biomass) were 

collected for 281 plants (4—22 plants per target site and site of 

1 origin). Finally, mycorrhizal colonisation was assessed for 191 
S" plants (68% of the surviving plants) because of the poor quality of 
J the roots in the remaining plants (roots that were too thick or too 
g short). 

a 

Data analysis 

g To analyse the data in the reciprocal transplant experiment, we 

used an analysis of variance for normally distributed variables 
§ (mycorrhizal colonisation, the number of leaves and aboveground 

Q, biomass) and an analysis of deviance for binomial variables 

(survival). We log-transformed both the number of leaves and the 
3 aboveground biomass to fit the assumptions of normality. We 

E distinguished the target region and site (where the plants were 

^ planted experimentally) and the region and site of origin (where 

5 the plants had come from) of the transplanted plants. We tested 

^ the target region against target site and target site against block 

S (five blocks, represented by rows in each site; Table 3). Similarly, 

^ we tested region of origin against site of origin and site of origin 

o against the target site-by-site of origin interaction (Table 3). Then, 

= we tested the interaction between target region and region of 

^ origin against the interaction between target site and site of origin, 

^ which was tested against residuals (Table 3). A significant 

"m interaction between target region and region of origin would 

0 indicate local adaptation at the regional scale. Similarly, a 
^ significant target site-by-site of origin interaction would indicate 
^- local adaptation (or maladaptation) at the smaller scale. We also 
^ tested for local versus foreign contrast, as was performed in [21], 
" but it was not significant and therefore we do not show the results. 
-° Where the interaction between the target region and region of 
_o origin was significant, we tested for the effect of region of origin in 
c each target region separately. We followed the statistical design of 
■g the study and tested the region of origin against the site of origin 
3 and the site of origin against residuals. 

.E To test for the relationship between mycorrhizal colonisation 

% and plant growth, we used a general linear model (GLM). We used 

c the number of leaves and aboveground biomass (both log scale) as 

g" dependent variables and the percentage of mycorrhizal colonisa- 

c tion as an independent variable. Because we expect that 

1 o mycorrhizal colonisation of the plants will differ among sites and 
^ regions, we used target region, target site, region of origin and site 

§ m of origin as covariates in the analyses. All tests were performed in 

the statistical program S-plus 6.2 (Insightful Corp., Seattle, 
Washington, U.S.A.). 



Q. (U 



g 5 



Results 

Both target site and block within the target sites significandy 
affected mycorrhizal colonisation, survival percentage, the number 
of leaves and aboveground biomass, although the target region did 
not (Table 3). The region of origin significantly affected 
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Table 3. Summary of the analysis of variance (F values) and deviance (Quasl-F values) for tlie effects of target region (Rt), target 
site (St; nested within target region), blocl< (nested within target site), region of origin (Ro), site of origin (So; nested within region 
of origin) and their interactions on survival percentage (Survival), the logarithm of the number of leaves (Leaves), the logarithm of 
the aboveground biomass (Biomass) and percentage of mycorrhizal colonisation (Colonisation). 





DF 




Survival 


Leaves 


Biomass 


Colonisation 








Quasi-F 


F value 


F value 


F value 


Rt 


1 


St 


0.73 


0.41 


0.28 


1.22 


St 


2 


Block 


6.90** 


12.15*** 


3.83* 


6.91** 


Block 


16 


Residuals 


2 59*** 


2.03* 


4.98*** 


4.19*** 


Ro 


1 


So 


1.46 


5.14+ 


26.08** 


49.86** 


So 


4 


St X So 


1.45 


2.22 


0.26 


0.56 


Rt X Ro 


1 


St X So 


0.43 


4.15+ 


0.40 


0.31 


St X So 


14 


Residuals 


1.35 


0.60 


0.80 


1.66+ 


Residual df 






507 


241 


241 


151 



Error indicates the error term used for each source of variation. 
***p<0.001; **p<0.01; *p<0.05; +p<0.1. 
doi:l 0.1 371 /journal.pone.0093967.t003 

mycorrhizal colonisation and the aboveground biomass, and it also 
affected the number of leaves with marginal significance (Table 3). 
Specifically, plants originating from region S had higher mycor- 
rhizal colonisation, higher aboveground biomass and more leaves 
than plants originating from region K (Table 2; Fig. 1). 

There was a marginally significant difference in mycorrhizal 
colonisation between the sites of origin within target sites (St x So 
interaction; Table 3). The number of leaves differed slightly 
between regions of origin within target regions (Rt x Ro 
interaction; Table 3). In region S, plants from the local region 
had significantly more leaves than plants from the foreign region 
(Fi 4 = 1 1 .93, p = 0.026), while the number of leaves did not differ 
significantly between the regions of origin in region K (Fi 4 = 0.72, 
p = 0.44; Fig. lb). 

After adjusting for the effects of target site and site of origin, we 
found a significant positive relationship between mycorrhizal 
colonisation and aboveground biomass (GLM; log scale; r = 0.16, 
p = 0.013, df error = 169) but not the number of leaves (GLM; log 
scale; r = 0.03, p = 0.703, df error = 181). 

Discussion 

Mycorrhizal colonisation 

Through a reciprocal transplant experiment in the field, we 
demonstrated that mycorrhizal colonisation has a genetic basis 
and varies among regions of origin. In both regions, plants from 
region S had higher mycorrhizal colonisation than plants from 
region K. This finding agrees with previous greenhouse experi- 
ments using A. amellus [10,20]. Similarly, a genetic basis of 
mycorrhizal colonisation was found in a greenhouse experiment 
with Plantago lanceolata, in which the percentage of mycorrhizal 
colonisation in the experiment corresponded to the field coloni- 
sation [9]. In contrast, another study using Andropogon gerardii 
suggested that the plants can change their percentage of root 
colonisation depending on soil nutrient content [8] . Overall, this is 
the first study that showed that mycorrhizal colonisation has a 
genetic basis using a field experiment. Furthermore, we demon- 
strated that the predetermination to low and high percentage of 
root colonisation persisted after five years of growth in foreign 
conditions. 




K S 
Target region 



Figure 1. Differences in mycorrhizal colonisation and plant 
growth parameters between regions of origin and target 
regions, a) Mycorrhizal colonisation, b) The number of leaves, and c) 
Aboveground biomass. Region K is Czech Karst, region S is Ceske 
Stredohori. Values represent the means ± 1 SE. To allow ecological 
evaluation of the data, we present untransformed data for the number 
of leaves and aboveground biomass. Asterisk (*) denotes significant 
effect (p<0.05) of region of origin within a target region. 
doi:10.1371/journal.pone.0093967.g001 
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Local adaptation 

We found no evidence of local adaptation in terms of plant 
growth in this experiment at the population or regional level, i.e., 
no significant interaction between target site and site of origin or 
between target region and region of origin. There was only a 
marginally significant target region-by-region of origin interaction 
in the number of leaves. However, the number of leaves is not an 
appropriate measure of fitness in the absence of an effect on plant 
biomass. In long-lived perennials, belowground biomass may also 
be an important indicator of plant fitness, as is aboveground 
biomass. However, it is not possible to completely excavate 
belowground biomass in the field. Therefore, we could not 
estimate this important parameter in our study. Nevertheless, 
some studies have shown a strong correlation between above- and 
belowground biomass in garden experiments (e-g. [7,10]). 

Generally, plants from region S outperformed plants from 
region K in both target regions. In a previous study that compared 
the performance of A. amellus in the same reciprocal transplant 
experiments after two years, we found some evidence of local 
adaptation at the population level [21]. However, local adaptation 
was found only in the seedling stage and differed among 
pojjulations [21]. As in the current study, no local adaptation 
was found in the transplant experiment for juvenile plants after 
two years [21]. 

The absence of local adaptation has also been reported in other 
plant species, such as Carlina vulgaris [32], Crotalaria pallida [33], 
Stipa capillata [34] and Helleborus foetidus [35]. Each of these studies 
provided different explanations for the lack of local adaptation. 
The absence of local adaptation in Carlina vulgaris at the local 
spatial scale, despite strong genetic difiFerentiation between 
populations, was explained by potential adaptation to other 
factors (such as pathogens or mutuaUstic organisms) that were not 
included in the study [32]. The lack of regional adaptation in the 
legume plant Crotalaria pallida was attributed to the unstable 
population structure of its herbivore [33]. Furthermore, the 
absence of local adaptation of Stipa capillata to their soil biota from 
Europe and Asia could be due to the similar composition of 
rhizosphere biota across populations and regions, sufficient gene 
flow, selection in favour of plasticity, or a functional redundancy 
among different soil biota [34] . Moreover, the lack of evidence for 
local adaptation in Helleborus foetidus could be caused by a 
congruency in selective pressures exerted by different soil 
environments on seedling emergence and survival [35]. 

In our study, plant aboveground biomass and mycorrhizal 
colonisation exhibited corresponding differences between the two 
target regions and the regions of origin. Because all of the studied 
sites are rather nutrient-poor, arbuscular mycorrhiza should have 
beneficial effects on the growth of Aster amellus according to the 
coadaptation model [12,14]. A previous experiment, where the 
percentage of mycorrhizal colonisation was manipulated, demon- 
strated the positive effects of mycorrhizal colonisation on the 
growth of A. amellus (Pankova et al., unpubl.), a finding that 
supports the coadaptation model. In the present study, we did not 
find any significant relationship between mycorrhizal colonisation 
and the number of leaves, and the positive relationship between 
mycorrhizal colonisation and aboveground biomass was quite 
weak (r = 0.16). Nevertheless, we tested for the relationship after 
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adjusting for the effects of target site and site of origin. Therefore, 
the weak relationship could be explained by a low variation in 
mycorrhizal colonisation within natural populations [20]. As 
mycorrhizal colonisation differs significantly between the two 
regions of origin, the main effect of mycorrhizal colonisation on 
the growth of A. amellus was pronounced at the regional scale. 

The genetic differences in mycorrhizal colonisation between 
plants from the two regions under study indicate that the 
percentage of mycorrhizal colonisation in the roots is an adaptive 
trait. Particularly, mycorrhizal colonisation was higher in the 
plants from region with lower soil nutrient content (region S). It is 
likely that higher mycorrhizal colonisation is the reason why the 
plants from region S outperformed the plants from region K in 
both target regions. This implies that the plants from region K are 
maladapted via their low inherent mycorrhizal colonisation levels. 

Conclusions 

Our study showed significant differences in mycorrhizal 
colonisation between regions of origin after five years in the field, 
indicating that the percentage of mycorrhizal colonisation of the 
roots has a genetic basis. This study confirms our previous findings 
from pot experiments using the same plant species A. amellus 
[10,20]. Overall, we did not find any evidence of local adaptation 
in terms of aboveground biomass in A. amellus in this study. 
Instead, plants from region S outperformed plants from region K 
in both target regions. Similarly, plants from region S showed 
higher mycorrhizal colonisation in all cases, which was driven by 
the lower nutrient content in the soil from this region. Thus, the 
pattern of plant aboveground biomass corresponded with the 
differences in mycorrhizal colonisation between the two regions. 
Higher mycorrhizal colonisation in the plants from region with 
lower soil nutrient content (region S) in both target regions 
indicates that mycorrhizal colonisation is an adaptive trait. 
However, lower aboveground biomass in the plants with lower 
mycorrhizal colonisation suggests that the plants from region K 
are in fact maladapted by their low inherent mycorrhizal 
colonization levels. We conclude that including mycorrhizal 
symbiosis in local adaptation studies may increase our under- 
standing of the mechanisms by which plants are adapting to their 
environment. 
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